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Abstract 

The acoustic peak in the CMB power spectrum is sensitive to causal 
processes and cosmological parameters in the early universe up to the 
time of last scattering. We provide limits on correlated spatial varia- 
tions of the peak height and peak position and interpret these as con- 
straints on the spatial variation of the cosmological parameters (baryon 
density, cold dark matter density and cosmological constant as well as 
the amplitude and tilt of the original fluctuations). We utilize re- 
cent work of Hansen, Banday and Gorski (HBG) who have studied 
the spatial isotropy of the power spectrum as measured by WMAP by 
performing the power spectrum analysis on smaller patches of the sky. 
We find that there is no statistically significant correlated asymmetry 
of the peak. HBG have also provided preliminary indications of a pre- 
ferred direction in the lower angular momentum range( t ~ 2-40) and 
we show how possible explanations of this asymmetry are severely con- 
strained by the data on the acoustic peak. Finally we show a possible 
non-gaussian feature in the data, associated with a difference in the 
northern and southern galactic hemispheres. 



1 Introduction 



One of the foundations of our description of the universe is the principle of 
isotropy and homogeneity - the statement that the universe is the same in 
all directions and locations. The isotropy of the universe can be explained 
by the theory of inflation, in which all of the observed universe descends 
from a single causally connected patch of the early universe. However, in- 
flation also contains the ingredients for generating anisotropies on various 
scales. This is manifest most importantly in the temperature anisotropies 
of the cosmic microwave background (CMB). The temperature anisotropies 
are generated by the quantum fluctuations of a scalar field during the in- 
flationary epoch. Other scalar fields could also have quantum fluctuations 
and these could in principle influence the distribution of matter in different 
portions of the universe. It is possible that these could leave a residual "tilt" 
to the universepp. 

Recently Hansen, Banday and Gorski (HBG) [2111311] have explored novel 
tests of isotropy by studying the WMAP 5 power spectrum extracted from 
the CMB using patches in different directions of the sky. This is an impor- 
tant test because the CMB provides the largest lever arm for the study of 
spatial asymmetries. In fact, in the lower angular range of i ~2-40 HBG 
find preliminary indications (at about the 2-3 sigma level) for non-zero dif- 
ferences in the power spectrum in hemispheres oriented along galactic co- 
latitude and longitude (80°, 57°), close to the ecliptic pole. This directional 
asymmetry seems to be connected with the orientation of the quadrupole 
and octopole p]. If this signal proves to be real, it could be an important in- 
dication of a fundamental tilt in the properties of the universe. In addition, 
a lack of isotropy has been seen in tests for non-gaussian behavior. Eriksen 
et al [7] and Park [S| were the first to find a difference between the north- 
ern and southern galactic hemispheres in tests for non-gaussianity, and this 
hemisphere difference has subsequently been confirmed by other measures 

Hams muni- 

HBG also studied the height and the location of the first acoustic peak 
on different patches of the sky, finding no obvious anisotropy within the 
expected cosmic variance. We use the HBG data to study the properties of 
the acoustic peak in more detail. The region of the acoustic peak is special 
because it is sensitive to causal processes in the early universe up to the time 
of last scattering. These processes depend on the cosmological parameters 
and any gradient of the parameters would yield a specific pattern between 
the peak height and peak position, i.e. both the height and the position 
would change in a correlated way. These correlations could reveal a sig- 
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nal of anisotropy that is not visible when looking at the height or position 
individually. We explore the HBG data for the possible presence of these 
correlated variations. Interpreting the results as variations in the cosmolog- 
ical parameters, we provide upper bounds on their possible variation. When 
we allow for an arbitrary correlated variation, we will find a weak correla- 
tion to a directions on the sky, but we will argue that it is not statistically 
significant. We also use the properties of the acoustic peak as a probe of 
the possible physics that could produce the low I signal found by HBG. The 
isotropy of the peak is a strong constraint on the physics that could produce 
the low t signal. If one studies only the variation of a single cosmologi- 
cal parameter, only a difference of the spectral tilt could explain the low-£ 
anomaly while being consistent with the properties of the acoustic peak. 
Finally we describe an unusual behavior in the distribution of the proper- 
ties of the acoustic peak, associated with a difference between the northern 
and southern galactic hemispheres. This may possibly be related to the 
north-south difference which has been found in tests for non-gaussianity. 



2 Characteristic patterns of individual variations 

The acoustic peak is categorized by its location in i and by its height. How- 
ever, these are not fully independent variables in terms of their dependence 
on the cosmological parameters. If one were modify the baryon density, for 
example, there would be a correlated variation in both the height and loca- 
tion of the peak. In this section we describe the nature of these correlations. 

Consider the variation of the cosmological parameters. In all cases, we 
adopt as the central value that found by the WMAP team. In addition, we 
define a parameter 

„. _ „WMAP 

e -h n m 

tj WMAP v-v 

where pi is any of the cosmological parameters and the superscript WMAP 
denotes its value measured by WMAP. 

The dependence of the power spectrum on all the cosmological parame- 
ters can be calculated using CMBFAST 13j. Here we focus on the properties 
of the acoustic peak. It is important to use the version of CMBFAST that 
does not implement the COBE normalization. Otherwise, changes in the 
power at low values of t would have an extra impact on the amplitude of the 
acoustic peak by rescaling the overall power spectrum. This is especially rel- 
evant for variations of the spectral tilt and the cosmological constant, which 
appreciably modify the power at lower values of I. Taking small variations of 
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the cosmological parameters with modify the location of the acoustic peak, 
£ p , and the height of the peak, h p , from the values found by WMAP, called 
Iq and ho. In each case we keep Otot = 1 as we vary the parameter. For 
small variations the changes will be linear in e^, so that we parameterize the 
results by 

t P = io + ai€i 

hp = h + bi€i (2) 

Thes can be combined to describe the correlation characteristic of each pa- 
rameter, using 

tp -to = n{h p - ho) (3) 

where 

ri = f (4) 

The value of the different coefficients are given in Table 1. The correlated 
variations in the (h, I) plane are shown in Figure 1 for the different param- 
eters. 



parameter 


Baryon 


CDM 


A 


amplitude 


tilt 


a 


10 





-22 


-16 


37 


b 


1875 


-2381 





5600 


-5491 



Table 1: Parameters describing the variation of the height and location 
of the acoustic peak when one changes the energy density of baryons, the 
energy density of cold dark matter, the value of the cosmological constant 
or the value of the spectral tilt. 



3 Measures of Isotropy 

Hansen et. al. (HBG) 2 have studied the WMAP power spectrum on 34 
patches of the sky, each with an opening angle of 19°. These cover most of 
the sky, avoiding the galactic plane where foreground contaminations could 
distort the signal. The HBG data on the location and height of the acoustic 
peak is displayed in Fig 1. The correlation with position on the sky is given 
in Table 2. The expected spread in these values has been estimated by HBG 
through the variance found in 512 simulated skies with gaussian distributed 
temperature fluctuations using the WMAP power spectrum. This yields a 
variance in the i direction of oi = 6.5 and in the height of = 300. 
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Figure 1: The lines show the correlated variation of the height and location 
in £ of the acoustic peak as one varies each of the indicated quantities. Also 
shown are the HBG data for the height and location of the peak for 34 
patches of the sky. 



A comparison of the average values of the position and height of the 
peak over opposite hemispheres is subject to large fluctuations as individual 
patches near the borders of the hemispheres are counted first in one hemi- 
sphere and then in the other as the orientation of the hemisphere changes 
slightly in direction. Therefore we seek a measure comparing directions that 
is not as sensitive to the border between the hemispheres. To this end we 
consider a possible dipole asymmetry. Theoretically a dipole is often the 
leading multipole for a spatial asymmetry [H [T^J [TH1 HH1 1^1 HH1 1201 EI] . A 
dipole tilt in the Universe can also be hidden beneath the overall dipole seen 
in the CMB temperature, which is usually attributed to our motion with 
respect to the CMB, but which can also have other components [T| IT5] . 

Consider a reference direction given by polar angles 8, <j) and let 7 be 
the angle between this reference direction and a point in the sky. For a 
given parameter pi, a correlated dipole variation centered on this reference 
direction has the form 

£(9,4>) = £0 + o-if] cos 7 

h(0,(f>) = h + 6^ cos 7 (5) 
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where a^, bi are the constants determined in the previous section and r] is 
the amplitude of a possible variation. 

For each of the cosmological parameters, we allow the reference direction 
to vary and search for a direction that gives the maximum possible dipole 
asymmetry. 1 While in each case a preferred direction can be found, the 
resulting asymmetries are not statistically significant, as we discuss below. 
The amplitude and standard deviation for each variation is shown in Table 
3. For each of the parameters, we can convert our lack of a signal into a 
bound. For the direction that maximizes the asymmetry we take the 2a 
upper bound on the amplitude 77. The results are given in Table 3. 

We also searched for a dipole anisotropy with an arbitrary correlation 
between the height and the peak location. We do this by allowing the relative 
size of a,i and bi to be arbitrary. The result was an anisotropy centered on 
the direction (9 = 107° ± 34°, = 30° ± 42°) with a small amplitude that 
is nominally 1.6 standard deviations from zero and a correlation 

r = ^ = -0.3 ±2.0 (6) 

The ratio r has a large uncertainty but favors a signal which is mostly in 
the overall height of the amplitude. Numerically, the amplitude corresponds 
to bi] = 148 ± 93, which is a signal of less than 3% and a nominal statis- 
tical likelihood of 1.6 standard deviations. Although this likelihood is of 
marginal significance, the true significance is weaker because we have al- 
lowed an arbitrary correlation by accepting any ratio r. To demonstrate 
this, we performed a similar search on 100 sets of simulated data, generated 
randomly with a gaussian weight and the same standard deviations given by 
HBG. If we allow the signal to have an arbitrary value or r and an arbitrary 
direction, we find that 85% of the simulated data sets also have a non-zero 
amplitude of at least 1.6 standard deviations. We conclude that the signal 
is not statistically significant. 

4 Constraints on the possible asymmetry at low £ 

In the previous section, we studied the constraints imposed by the proper- 
ties of the acoustic peak. However, some of the parameters can also modify 
the power spectrum at lower values of I. Indeed a study at lower values is 
especially interesting because HBG have identified a potentially significant 
asymmetric signal for the values t ~ 2 — 40, with the hemisphere oriented in 

1 These studies were performed with MINUIT. 
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the directions (80°, 57°) having about 20% more power than in the opposite 
hemisphere. Subsequent work by Hansen et.al. (HBBG)[Hj has explored the 
possibility that this signal could be due to correlated changes in the com- 
bination of the optical depth and the initial amplitude of fluctuations, or 
alternatively in the combination of the amplitude and tilt of the initial fluc- 
tuations. In this section, we explore how the constraints from the acoustic 
peak impact on possible explanations of the asymmetry found by HBG. 

Of the cosmological parameters that we study, the ones that influence the 
power spectrum at lower values of £ are the overall amplitude, the spectral 
tilt and the cosmological constant. The latter is by far the easiest to rule 
out as the source of the HBG asymmetry. A variation in the cosmological 
constant would lead to changes at low values of £ through the late time 
integrated Sachs- Wolfe effect and also lead to changes in the acoustic peak 2 . 
However, in order to produce a 20% change in the power spectrum at low 
£, one would require over 100% variation of the cosmological constant. Our 
constraints given above firmly rule this out. 

The simplest variation that would explain the HBG signal is that of the 
overall amplitude, requiring a 20% larger amplitude in one hemisphere than 
the other. Since the whole power spectrum is proportional to this amplitude, 
the same effect would show up at the acoustic peak. This is not seen, as is 
clear from our overall constraint of the previous section. Moreover, when 
we study the amplitude variation centered on the HBG direction, we find 
an even tighter constraint, with a best fit amplitude of rj = 0.02 ± 0.02 and 
a 2cr bound of r] < 0.06. 

Finally a variation in the tilt of the CMB spectrum would also influence 
the relative height of both the low and high £ regions. Here we are guided by 
the final result of HBBG in the section where they consider the possibility 
of correlated variations of both the tilt and the amplitude. Their central 
values are consistent with the same amplitude but a tilt that is different 
in the two hemispheres by 8%, a 2-3a effect. We test this hypothesis by 
considering the properties of the acoustic peak under such variation. The 
spectral tilt has the property that it influences the power spectrum more 
strongly at lower values of £ than it does at the acoustic peak. Our studies 
with CMBFAST reveal that a 20 ± 8% variation in the power at low values 
can be fit with a variation of the spectral index of 0.07 ± 0.03, and that 
this would be expected to produce a variation in the height of the acoustic 
peak of 7 ± 3%, with a corresponding change in the peak location. We 

2 The cosmological constant does not influence the causal physics up to the time of last 
scatter, but it influences the angular scale of the acoustic peak. 
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can search for such a variation, finding a difference of 3 ± 2%, i.e. a small 
but not statistically compelling signal. However, within the error bars, this 
observation is compatible with the expected signal. 

We also can consider the combined variation of both the amplitude and 
the spectral tilt, using the information at lower values of i as well as infor- 
mation at the peak. Performing a simultaneous fit to both regions, we find 
consistency for 

n N -n s = 0.083 ±0.041 

4^ = 1.05 ±0.025 (7) 

This result is consistent with that of HBBG within the error bars. 



5 A hint of non-gaussianity 

We have found no significant signal of a correlated variation of the height 
and location of the acoustic peak. However the data does appear to contain 
a more unusual asymmetry, which we discus in this section. The northern 
and southern galactic hemispheres appear to have a different variance in the 
spread of the heights and locations of the peak, with the northern hemisphere 
having more variation. This is visible to the eye in Fig. 2. To probe this 
effect we define a variable 

5i =J{^—k)2 + { h^)2 (8) 

The distribution of the values of 5 for the northern and southern hemispheres 
is shown in Fig. 3. It is easy to see visually that there is a significant 
difference between the hemispheres, with the northern hemisphere appearing 
to be anomalous. We refer to this as a non-gaussian feature since at would 
imply that there is not a uniform gaussian distribution in the properties of 
the peak. 

We search for the direction that has the greatest non-gaussianity by 
identifying the direction that maximizes the asymmetry 

/(M) = J>cos 7 i ( 9 ) 

i 

This search identifies 9 ~ n as the direction that maximizes the difference, 
i.e. there is an asymmetry between the northern and southern galactic 
hemispheres. The net difference corresponds to f{9 = 7r, </>) = —7.5. 
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Figure 2: The data points on the acoustic peak are here separated into 
those of the northern (open circles) and southern (filled circles) galactic 
hemispheres. One can see visually that the data for the northern hemisphere 
is more widely scattered than that of the southern hemisphere. 



We study the significance of this signal by simulating 256 data sets chosen 
randomly with a gaussian weight and the standard deviations given by HBG. 
In each case we search for the direction that maximizes the asymmetry 
f(9, (j)) and find how large the asymmetry is when centered on this direction. 
The simulations reveal a median of f(9, <fi) = —3.3 and a standard deviation 
of <jf = 1.3. Consistently, we find that only one of the 256 simulated data 
sets has an asymmetry as large as is seen in the real data. We see that our 
signal is unlikely to be purely a statistical fluctuation, with an exclusion at 
the 3a level. 

If this signal is real physics, it appears as an extra random ingredient 
in the temperature fluctuations - dominantly in the northern hemisphere. 
Since it involves regions that were out of causal contact at the time of last 
scatter, it would either correspond to local foreground effects in the present 
universe or to physics very early in the period of inflation when these regions 
were causally connected. While the signal that we have identified appears to 
be statistically significant, it is not clear how it is related to other measures 
of non-gaussian behavior. The use of other measures has produced ample 
evidence for non-gaussianity in the WMAP datapH El El EB EEU H21 E3 E31 
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Figure 3: The distribution of the parameter 5 in the northern (light) and 
southern hemispheres (dark). The distribution for the southern hemisphere 
is consistent with the expected gaussian distribution within statistics. As 
described in the text, the distribution for the northern hemisphere appears 
more dispersed than expected. 



1241 1251 1261 1271 1281 I29j . Moreover, as mentioned in the introduction, for 
several of these signals there seems to be a correlation with the northern 
and southern hemispheres [71151 151 1 1 Ul ITT1 1 1 2 j . For example, using a measure 
of the frequency of hills, lakes or saddles in the temperature distribution, |12j 
found that the northern hemisphere is non-gaussian at the 2-3 sigma level, 
while no significant non-gaussianity was found in the southern hemisphere. 
Our signal shares this directionality, even including the possibility that it is 
the northern hemisphere that is anomalous. 

6 Conclusion 

We have studied the properties of the acoustic peak in the CMB power 
spectrum and have provided measures of isotropy for the types of correlated 
variations of the height and location of the peak. Within the statistical 
uncertainty, no significant asymmetric signal was found. 

The properties of the acoustic peak also constrain the possible explana- 
tions of the tentative signal of an asymmetry at low i found by HBG. If 



9 



one attempts to explain this by the variation of a single parameter from the 
set considered in this paper, only the variation of the spectral tilt is able 
to accommodate the constraints of the acoustic peak. This is in agreement 
with the overall fit performed by HBBGjS]. 

We have also identified a curious difference between the data in the 
northern and southern hemispheres, with the northern hemisphere being 
more widely scattered than expected from the cosmic variance. 
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co — latitude 


longitude 


I location 


amplitude 








223.3 


5551 


35 





217.0 


5235 


35 


60 


215.8 


5515 


35 


120 


218.1 


5909 


35 


180 


207.5 


6211 


35 


240 


231.8 


5231 


35 


300 


232.7 


5758 


65 





212.3 


5957 


65 


36 


217.9 


4640 


65 


72 


230.9 


5410 


65 


108 


224.5 


5099 


65 


144 


223.5 


5506 


65 


180 


221.5 


5888 


65 


216 


213.1 


5506 


65 


252 


217.9 


5986 


65 


288 


222.3 


5927 


65 


324 


224.6 


5335 


115 





219.5 


5633 


115 


36 


222.3 


5415 


115 


72 


221.3 


5510 


115 


108 


219.5 


5615 


115 


144 


224.8 


5804 


115 


180 


216.3 


5693 


115 


216 


230.5 


5627 


115 


252 


220.8 


5113 


115 


288 


225.9 


5395 


115 


324 


219.4 


5759 


145 





218.8 


5315 


145 


60 


219.3 


5683 


145 


120 


221.9 


5792 


145 


180 


216.6 


5302 


145 


240 


217.1 


5362 


145 


300 


221.4 


5636 


180 





218.7 


5513 



Table 2: HBG data on the height and location of the acoustic peak on 
patches centered on different locations on the sky. We thank the authors of 
[2] for permission to include this data. 



13 



parameter 


Baryon 


CDM 


A 


amplitude 


tilt 


j] fit 


0.076 ± 0.048 


0.004 ±0.11 


0.055 ± 0.035 


0.027 ±0.017 


0.03 ±0.02 


2a bound 


0.17 


0.28 


0.125 


0.061 


0.057 



Table 3: Isotropy bounds when interpreted as variations of the underlying 
parameters. The first row of results gives the fit value of the fractional 
variation of each parameter, while the second row reports the 2 a upper 
limit on the possible variation. 
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